Finite size effects in the dynamics of opinion formation 
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For some models of relevance in the social sciences we review some examples in which system 
size plays an important role in the final outcome of the dynamics. We discuss the conditions under 
which changes of behavior can appear only when the number of agents in the model takes a finite 
value. Those changes of behavior can be related to the apparent phase transitions that appear in 
some physical models. We show examples in the Galam's model of opinion transmission and the 
Axelrod's model of culture formation stressing the role that the network of interactions has on the 
main results of both models. Finally, we present the phenomenon of system-size stochastic resonance 
by which a forcing signal (identified as an advertising agent) is optimally amplified by a population 
of the right (intermediate) size. Our work stresses the role that the system size has in the dynamics 
of social systems and the inappropriateness of taking the thermodynamic limit for these systems. 



I. INTRODUCTION 

In statistical physics we are used to taking routinely 
the thermodynamic limit in which the number of con- 
stituents (or, more precisely, the number of degrees of 
freedom) N tends to infinity 1]. This is necessary in 
many aspects, such as determining the validity of the 
statistical approach, the equivalence between ensembles, 
the existence of phase transitions with symmetry break- 
ing and discontinuities in the order parameters or their 
derivatives, etc. When we apply the methods of sta- 
tistical physics to explain the properties of macroscopic 
matter, it is clear that the number N is always finite, but 
very large. A good estimative of its order of magnitude is 
the Avogadro number N = 6.023 x 10 23 . Consequently, 
the nowadays widely used computer simulations of phys- 
ical systems always struggle to get to larger and larger 
systems with the continuous increasing demand in com- 
puter resources. In many cases as, for instance, when 
trying to predict the thermodynamical properties near 
a critical point, it is essential to use the well developed 
techniques of finite size scaling 2] in order to extrapolate 
the results to A~ — > oo. 

When applying the same tools of statistical physics to 
problems of interest in social sciences 0. [J. 13 ■ we have 
to take into account that the number of individuals or 
agents considered can never be that large. In most cases, 
realistic values of N range in the hundreds or thousands, 
reaching at most a few million. The thermodynamic limit 
might no be justified in this the results in that 

limit can vary with respect to those of finite-size systems. 
Furthermore, new and interesting phenomena can appear 
depending on the number of individuals or agents con- 
sidered. In this paper we will review some finite-size re- 
lated effects that have been recently found in some mod- 
els of interest in social systems and that exemplify cases 
in which the thermodynamic limit does not capture the 
essentials of the system's behavior. 

The outline of the paper is as follows: in the follow- 
ing section III Al we will explain the concept of apparent 
phase transitions and pseudo-critical points by reviewing 



a physical example in which a clearly observable change 
of behavior disappears in the thermodynamic limit. We 
will then show in the following sections that apparent 
phase transitions also appear in two widely used models: 
the Galam model for opinion spreading (in section Til B|) 
and the Axelrod model of culture formation (in section 
III C|) : finally, section ITTT1 is devoted to the phenomenon of 
system-size stochastic resonance in a model for opinion 
formation that takes into account external influences. A 
brief summary and concluding remarks are presented in 
the final section ITVl 



II. APPARENT PHASE TRANSITIONS 
A. Pseudo-critical temperatures 

A true phase transition can only appear in the ther- 
modynamic limit 0. This limit is needed for the singu- 
larities inherent to the transition (for instance, the di- 
vergences in the critical point of a continuous, second 
order, phase transition or the discontinuities in a first or- 
der phase transition) to develop from the sum of analytic 
functions. However, in some cases, finite systems show 
a behavior that reminds us of a phase transition. An 
example which has been extensively studied is the sine- 
Gordon model for the growth of surfaces [|| Q • Nowa- 
days, it is well established that in two or more spatial di- 
mensions the model displays a Kosterlitz-Thouless-type 
phase transition from a (low temperature) smooth phase 
to a (high temperature) rough phase at a g iven critical 
value T c of the temperature H, & El El EL 

However, numerical simulations of this model seemed 
to imply that the same transition was also present in 
the one-dimensional case |l3j . See the evidence in figure 
n where it can be observed that the roughness of the 
one-dimensional surface is zero until a temperature T c ks 
1 and then increases suddenly. Despite this numerical 
evidence, the result was doubtful, since the theory clearly 
indicated that such a phase transition was not possible 
for a one-dimensional system 0, ITM . In the view of 
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FIG. 1: Roughness lu 2 /N of the one- dimensional sine-Gordon 
model as a function of the temperature 13] . Observe the ap- 
parent transition between a smooth phase (roughness equal 
to zero) to a rough phase at a pseudo-critical temperature 
T c ~ 1. Different symbols correspond to different system 
sizes, as indicated in the figure and the line is the result 
uj 2 /N = T/12 of the Edwards- Wilkinson model valid at high 
temperatures fl7|| . 



this negative prediction, one can revisit the simulation 
data and conclude that the deviation from smooth 
to rough surfaces occurs at a pseudo-critical temperature 
that depends on the system size as T C (N) ~ l/ln(iV). 
It is certainly true that in the thermodynamic limit the 
transition point switches to T c = 0. In other words, the 
system does not display a phase transition, as demanded 
by the rigourous theorems. However, it is also true that 
for every finite system, a well defined transition point 
(which is not a true critical one) can be identified as 
separating the smooth from the rough behavior. Since 
the dynamics at low temperature is dominated by the 
existence of many metastable states with large energy 
barriers jS] , a lot of numerical effort is needed in order to 
proof beyond any doubts that the transition point tends 
to zero for increasing system size|l3j. 



B. Galam's model 

There are several models that a naly ze the dynamics of 
opinion formation in a socie ty fi g, jlfl. l20 | , l2lL p2l [2a, 

Elj . The model introduced by S. Galam |2(J,[21J aims to 
identify mechanisms by which an opinion, being initially 
minority, can become majority in the long run. There 
have been many recent examples of such change of opin- 
ion: the rejection to the Maastricht treaty in Ireland, the 
negative vote in France for the European constitution, 
the switch in the belief of the authorship of the March 
11th 2004 terrorists attacks in Madrid, the rumors con- 
cerning some September 11th opinions in France plj . etc. 
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FIG. 2: The inset shows p as a function of p and N for 
Galam's model in which groups form independently of any 
newtork structure. The size distribution of discussion groups 
is uniformly distributed in [1,M]. Open symbols correspond 
to M = 4 (p c = 0.2133), while filled ones to M = 5 
(p c = 0.3468). Different symbols correspond to different 
system sizes between N — 10 3 and N = 10 6 (the tran- 
sition sharpens with increasing system size). The main 
plot shows the collapse according to the scaling function 
p(p,N) — (1 + erf(a;/1.17))/2 as a function of the scaling 
variable x = (p — p^N 1 ^ 2 . 



The basic ingredient of the model is that the opinion a 
social group has on a topic changes when the individu- 
als meet to discuss in small groups. To simplify matters, 
the model considers that the opinion of an individual 
can only take two possible values: in favor or against. 
After discussion, the majority opinion is adopted by all 
the members of the discussion group. In case of a tie, 
i.e. same number of individuals against and in favor, one 
of the two possible decisions is always favored and hence 
adopted. This favored decision usually responds to some 
sort of social inertia against change. For instance, when 
discussing about the reforms in the European constitu- 
tion, the opinion favored by the French people was that 
of voting against its approval because it was unclear the 
advantages that the new constitution would bring to the 
French society. The model has been recently modified in 
order to acount for the diffusion of the individuals [3l| 
and the effects of contrarians (people holding a different 
opinion of that of the majority) [42ll43|. 

If one assumes that the discussion groups are formed 
by taking at random individuals from the whole soci- 
ety, it is easy to establish a recurrence relation for the 
fraction of people P+(i) holding the favored opinion at 
time t. This assumption by which any two individuals 
have the same chance to meet in a discussion group be- 
longs to the category of mean-Held approaches, widely 
used as a first ap prox imation in many problems of in- 
teracting particles[44j. The recurrence relation adopts 
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the general form: P + (t + 1) = F[P+(i)], where F is some 
polynomial function that depends on the size distribution 
of the discussion groups [21| ■ The analysis of this recur- 
rence relation shows that there are three fixed points: the 
trivial ones (which are stable) P + = and P + = 1 , and a 
non-trivial faith point P+ = p c (an unstable fixed point). 
Let us denote by p the initial proportion of supporters 
p = P+(0) of the favored opinion. If this proportion is 
greater than the faith-point value, p > p c , then the re- 
currence relation tends to the fixed point P + = 1 and 
the whole population adopts the favored opinion. Oth- 
erwise, if P+{0) < p c the fraction of people supporting 
the favored opinion tends to as the recurrence rela- 
tion proceeds in time. If we define an order parameter 
p = (P+(t — > oo)} as the average, over realizations of the 
discussion process, of the fraction of people holding the 
favored opinion after many discussion steps, the predic- 
tions of this simple analysis is that 

Since p c < 1/2 it is possible that a favored opinion that 
was held initially by a minority, p c < p < 1/2, is in the 
long run the opinion hold by the whole population. That 
would be the explanation, according to this model, of 
the fact that an initially minority opinion can become 
majority if the society is allowed to discuss about it for 
a long enough time. 

This analysis, based on the stability of the fixed points, 
is only valid in the N — > oo limit. In fact, for a small 
population, the condition p > p c docs not automatically 
ensure that the preferred opinion will be the favored one. 
There are finite size fluctuations such that the actual 
result of the repeated discussion process is not well es- 
tablished if the difference between the initial fraction of 
supporters and the faith point is of order \p— p c \ ~ N^ 1 / 2 . 
This rounding-off of the sharp transition can be summa- 
rized in the following finite-size scaling law: 

p(p,7V)-p((p-p c )iV 1 / 2 ) , (2) 

a result well confirmed by the computer simulations, as 
shown in figure [21 

It is possible to compute the time it takes the popula- 
tion to reach the final consensus. According to Galam's 
analysis, the recurrence relation reaches a de facto con- 
sensus after a finite, small, number of iterations. For a 
finite-size population, this can be quantified precisely |45| 
as the number of steps T required for P + (T) = 1 /N or 
P+(T) = 1 —1/N indicating that the vast majority of the 
population has adopted the same opinion. In figure [3] it 
is shown that this time scales as T ~ In AT. Numerical 
simulations, as can be seen in the same figure, also show 
that the time to reach consensus grows near the faith 
point; this slwowing down can be understood as initially 
the system escapes very slowly from the unstable fixed 
point p c . Therefore, the conclusion of this mean-field 
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FIG. 3: Time T to reach consensus as a function of N in the 
Galam's original model for different values of p and a max- 
imum group size of M = 5. Respectively from top to bot- 
tom (in the right end of the graph) the different curves show 
p = 0.20, 0.10, 0.05, 0.30, 0.60, 0.90. In the inset, the time T is 
plotted as a function of p for the same system parameter as in 
figure^ with the same symbols meaning than in that figure. 
Observe the slowing-down for p close to the critical points. 



analysis is that this model displays a true phase transi- 
tion between regimes dominated by one or the other of 
the two possible opinions. The transition is rounded-off 
for a finite system and finite-size scaling laws can be used 
to analyze the results in finite systems. 

More dramatic finite size-effects appear if we go beyond 
the mean-field theory and consider spatial effects within 
the so-called neighborhood models [45j. Contrary to the 
previous scheme in which an individual could join any 
discussion group, neighborhood models assume that dis- 
cussion groups are formed with some criterion of locality 
or geographical neighborhood. Some other spatial effects 
of this model have been discussed in references 0, ^7J • 
Here we consider that the N individuals are fixed at the 
sites of a regular lattice and the discussion cells are well 
defined closed, localized, rectangular regions. The exact 
dynamics is as follows: at each time step, an individual 
is chosen and a rectangle with length sides m x and m y 
is built around that individual. Both numbers, m x and 
m y , are chosen uniformly from an interval [1,M]. The 
m x x m y individuals on this region discuss about the se- 
lected topic and all reach a local consensus based upon 
the majority rule, with a bias (favored opinion) in case 
of a tie. The opinion of all individuals in the group is 
changed accordingly and time increases by an amount 
t —> t + m x x m,y/N. 

As in the original Galam's model, a global consensus is 
always reached in a finite time. The first feature that ap- 
pears in the numerical simulations of this neighborhood 
model is the smooth transition between the two different 
selected states as a function of the initial proportion p 
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FIG. 4: p as a function of p for different values of N in the 
neighborhood's Galam model. The maximum lateral size of 
the groups is M = 5. The system size varies from N — 225 
to N = 10 6 . The main plot shows the scaling law p(p, N) — 
p(pN a ) with a = 0.48. The inset shows the unsealed results, 
showing the transition point shifting toward p = as size 
increases. 



of supporters of the favored opinion. This is evident in 
the inset of figure where we plot the order parameter 
p as a function of p. This looks similar to the rounding- 
off induced by finite size of the original model. However, 
the new feature is that the transition point, defined as 
the fraction p for which p = 1/2, clearly decreases with 
system size N. In fact, a finite size scaling analysis shows 
that the data can be well fitted by the formula: 

p(p,N)=p( P N a ) (3) 

with an exponent a w 0.48 for the dynamical rules 
detailed above. The transition point scales hence as 
p c (N) ~ N~ a . An strict statistical mechanics analysis 
based on the thermodynamic limit N — > oo would con- 
clude that the actual critical point is p c = and hence 
the transition between the two regimes has disappeared! 
In our opinion, however, this would be a misleading con- 
clusion. There is indeed a well defined transition between 
the dominance of the two consensus opinions for a finite 
number N, of the order of magnitude to be encountered 
when applying this model to any situation of interest. 
In an infinite system, the conclusion is that the favored 
opinion eventually dominates regardless the (infinitesi- 
mally small) amount of initial supporters. In any finite 
system, though, it is possible to define a transition be- 
tween the two regimes at a value p c (N). Equivalcntly, 
for a fixed proportion of initial supporters p one could 
identify a value N c (p) separating the different consensus 
obtained whenever N > N c (p) or N < N c (p). Hence, 
the change in behavior is induced by the variation of the 
number of individuals, N. A result that would have been 
missed if we had taken the thermodynamic limit. 
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FIG. 5: Time to reach consensus as a function of TV in the 
neighborhood version of Galam's model, for p — 0.30 (circles) 
and p = 0.20 (squares). Two different scalig regimes for large 
and small N can be clearly seen. 



For the neighborhod model, it is possible also to com- 
pute the time T it takes to reach the final consensus. 
Again, the simulations indicate that this time increases 
near the above defined transition point p c (N), and this 
critical slowing down is a typical signature of the exis- 
tence of a phase transition. In this case, see figure 03 the 
characteristic time scales as T(N) ~ N@ with (3 ~ 0.6 for 
N < N c (p) and /3 - 0.1 for N > N c (p). The reason for 
these two regimes to appear is that, as discussed above, 
the final state to which the system approaches depends 
on system size: for finite p < 1/2 and small N the sys- 
tem will approach to P + = (as p = 0), but for a large 
enough system, the absorbing state will be P+ = 1, and 
thus the dependence of T with system size experiments 
a strong change. 



C. Axelrod's model 

In an important paper R. Axelrod has introduced 
a model for the generation and diversification of differ- 
ent cultures. The basic question addressed by Axelrod 
is the observation that differences between cultures do 
not disappear at all despite the fact that people tend 
to become more alike in their beliefs, attitudes and be- 
haviors when they interact. His proposal to explain this 
paradox is a model to explore mechanisms of competition 
between globalization (understood as a uniform culture) 
and polarization (understood as the persistence of cul- 
tural diversity). 

In the Axelrod model, culture is defined as a set of 
attributes subject to social influence. Each one of these 
attributes can take several traits. The traits change by 
social interaction between neighboring agents. The main 
premise is that the more similar an agent is to a neighbor, 
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the more likely the agent will adopt one of the neighbor's 
traits. In this way, the model incorporates the observa- 
tion that the communication is most effective between 
already similar people. The novelty this model brings 
into social modeling is that it takes into account the in- 
teraction between different cultural features. 

More specifically, Axelrod's model considers that each 
agent is characterized by a set of F cultural features, 
each one of them can adopt any of q different traits. The 
state of agent i is hence characterized by the F vari- 
ables (oil, . . . , CTjf), each one of them can take the val- 
ues <7jfe £ (0, . . . , q — 1). The number of possible traits 
per feature q can be considered as some sort of measure 
of the degree of heterogeneity a cultural feature might 
have. Agents are located in the sites of a lattice. The 
connectivity between the different lattice sites defines the 
neighbor structure of the social network. The dynamics 
is as follows: first, two neighbor agents i and j are ran- 
domly selected. One then computes their overlap £ij or 
number of common features, £,j = ^ fc 8 aik . ajk . With a 
probability £ij/F, the value of one of the not yet common 
features is transferred from one agent to the other, so in- 
creasing the overlap by one. This process is repeated by 
randomly selecting another couple of neighboring agents. 
Eventually, a frozen state is reached in which no possible 
evolution is possible. In such a frozen state, neighbor- 
ing sites have either an overlap equal to or to F. The 
relative size of the largest cluster of agents sharing all 
cultural features, S max /N is a measure of the cultural 
diversity. If S max = N it means that all agents share 
the same cultural features, a monocultural or globalized 
state. If S m ax — 1 it means that no neighboring agents 
share the same culture, a state that can be qualified as 
of cultural diversity or polarized. 

For fixed value of F it was shown by Castellano et at 
|35l l36j | that a genuine phase transition occurs between 
polarized and globalized states at a critical value, q c (F), 
of the parameter q. For q < q c {F) a monocultural state 
appears, whereas for q > q c (F) the society splits in sev- 
eral regions each one with a different culture. 

The analysis of Castellano et al. corresponds to the sit- 
uation in which fluctuations are neglected. The situation 
is dramatically altered when one considers the effect that 
cultural drift has on the dynamical evolution. Quoting 
Axelrod: "Perhaps the most interesting extension and at 
the same time, the most difficult one to analyze is cultural 
drift (modeled as spontaneous change in a trait)" 
We have modeled cultural drift by adding one step to 
the original Axelrod model: with probability r an agent 
changes randomly one of his cultural traits. For fixed 
values of F and q a change of behavior can be clearly ob- 
served when varying the rate r. A pseudo-critical value 
r c can be defined such that if r > r c the final state is 
multicultural, whereas for r < r c the cultural drift can 
not lead the system to develop many cultural states [73 . 
Note that we have used again the word pseudo-critical 
to denote the value r c . The reason, as before, is that r c 
depends on the system size N . The exact dependence 
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r N logN 

FIG. 6: For the Axelrod model denned in a regular two- 
dimensional lattice with nearest neighbors interactions, we 
plot the relative size of the largest cluster S max /N as a 
function of the probability of a random change in a cul- 
tural trait, r. Other parameters are F — 10, q = 100. In 
the inset, we plot the raw results for different system sizes 
N = 10 2 , 20 2 , 50 2 , 100 2 shown with circles, squares, diamonds 
and triangles, respectively. The main plot shows the rescaled 
results, showing a data collapse when the z-axis is rescaled as 
riVlniV. 



is related to the dimensionality of the underlying lattice. 
For a regular two-dimensional network of interactions, 
it is r c ~ 1/AlniV 38], see figure HH whereas in a one- 
dimensional lattice the dependence is r c ~ l/iV 2 [3jJ as 
seen in figure [7| In both cases, in the limit of infinite 
size iV — ► co, the value of r c tends to 0. In this limit, 
for any positive value of r it is r > r c and the system 
tends to a multicultural state independently of any other 
parameters. In the language of critical phenomena, the 
cultural drift is a relevant variable since any vanishingly 
small value alters the final state of equilibrium. However, 
in a finite system there is always a well defined value r c 
which separates the transition between polarization and 
globalization. Again, an analysis based upon the infinite 
size limit would have missed this interesting property of 
Axelrod's model. 

Another interesting modification of the Axelrod model 
considers complex networks of interactions between the 
agents. In order to mimic some of the conditions that 
a social network has, the connectivity is that of a small 
world network introduced by Watts and Strogatz [5(| : 
individuals are considered to be located in the sites of a 
regular network with nearest-neighbors interaction; with 
a probability p one of the nearest neighbor links is rewired 
and points to another, randomly selected individual, so 
establishing a long-distance connection. In this small- 
world network, the phase transition observed by Castel- 
lano et al. is still present, although its exact location 
shifts to larger values of the heterogeneity parameter q 
as the rewiring parameter p increases. This shows that, 
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FIG. 7: Same as figure [S] in the case of a one-dimensional 
lattice. The necessary scaling factor for the noise intensity is 
now r N 2 . The symbols correspond to different system sizes 

N = 10 4 , 10 3 , 10 2 shown with circles, squares and diamonds, 10 3 , ,5 x 10 3 , 10 4 shown with circles, squares, diamonds and 
respectively. Filled symbols use q = 5 while open symbols triangles, respectively, and F — 10. 
correspond to q — 10. In all cases it sis F = 10. 



FIG. 8: Axelrod model in a Barabasi-Albert lattice. The 
necessary scaling factor for the noise intensity is now r TV ' 39 . 
The symbols correspond to different system sizes N = 10 3 ,2x 

13 



as expected, the enhanced connectivity of the small- world 
network favors the globalized state. A different scenario 
appears in the scale- free networks. Those networks are 
characterized by a power law tail in the distribution of 
the number of neighbors. The first example is that of 
the Barabasi-Albert networkjHJ which is grown using 
the mechanisms of preferential attachment of links. In 
the Barabasi-Albert network, the location of the (pseudo- 
critical value q c depends on system size as q c {N) oc A^~ a 
with a ~ 0.39. This is a consequence of the finite-size 
scaling law by which the order parameter S max /N turns 
out to be a function of the product qN a , as evidenced 
in figure [SJ With the thermodynamic limit in mind, the 
conclusion is that the transition betweeen globalized and 
polarized states disappears in this case. However, as dis- 
cussed before, this conclusion does not prevent the finite 
system from having a true, observable, change of behav- 
ior as the parameter q varies. It should be said that the 
transition is recovered, even in the thermodynamic limit, 
when one considers other types of more complex networks 
such as the structured-scale-free networks 39| ■ 



III. SYSTEM SIZE STOCHASTIC RESONANCE 

It was found more than two decades ago that stochastic 
terms (noise) can have a constructive role in non-linear 
dynamical systems, |53l l54|. In these pioneer works, it 
was shown that the response of a dynamical system to 
a weak signal can be improved in presence of the right 
amount of noise, but when the strength of the stochastic 
terms is too large or too small, such improvement disap- 
pears. This phenomenon receives the name of stochastic 



resonance. 

The original motivation of these works was to provide 
a framework in which the small variations in the radia- 
tion reaching the Earth from the Sun could be amplified 
to account for the big climatic changes observed in ge- 
ological ages. However, stochastic resonance j5f| is 
now known to hold in a large variety of systems subjected 
to many different kinds of forcing. The paradigmatic ex- 
ample is that of an stochastic dynamics in a bistable po- 
tential under the effect of a weak periodic signal. In the 
absence of the forcing signal, the stochastic dynamics is 
such that there are jumps between the two stable wells of 
the potential with a characteristic Kramers time, 7V|57|. 
The role of the forcing is to lower periodically the bar- 
rier to go from one well to the other in such a way that 
it helps the stochastic dynamics to overcome that bar- 
rier alternatively in one direction or another following 
the forcing. If the value of half the period T/2 coincides 
with the Kramers time, then the external signal is opti- 
mally followed by the stochastic dynamics. The match- 
ing Tk = T/2 can be achieved by varying the intensity 
of the noise, in such a way that an optimal response to 
the action of the periodic forcing is obtained for the right 
amount of noise. 

The very general conditions under which the phe- 
nomenon can appear has opened a very active field of 
research. Many other situations have been found where 
noise plays a constructive role in a dynamical system. 
As in the prototypical example, all that is needed is a 
stochastic dynamics acting upon a system with two pos- 
sible equilibrium states and a weak periodic forcing. In 
fact, those conditions can be further relaxed, an interest- 
ing extension being that of coherence resonance by which 
an excitable or chaotic [5^ system can optimize the 



7 



regularity of its firings under the correct amount of noise. 
Similar results have been referred to as stochastic co- 
herence in or stochastic resonance without exter- 
nal periodic force |6ll . Another extension considers 
the enhancement of the response to non-periodic forcings, 
mostly within the context of some biological applications 

m 

Concerning the dependence on system size relevant to 
this review, Pikovsky and collaborators showed that 
an extended system can achieve a resonance effect with 
respect to an external signal as a function of the size of 
the system in a phenomenon that has been named as sys- 
tem size stochastic resonance. The explanation is that, 
due to the law of large numbers, the intensity of the noise 
experienced by the macroscopic variable scales as A~ 1//2 , 
being N a measure of the system size. Therefore, the ef- 
fective noise intensity for the, average, macroscopic vari- 
able can be controlled by varying the system size until 
the resonance condition between the Kramers time and 
the period of the forcing is satisfied. This very interest- 
ing result has been recently extended to the system size 
coherence resonance phenomenon, with a similar expla- 
nation |65| . 

In the dynamics of social systems, stochastic resonance 
was found in very simple opinion models [661 |6jj . The 
basic ingredients are (i) the individuals of a society can 
have a binary opinion (in favor or against) a particu- 
lar topic; this opinion can change by social interaction in 
such a way that both opinions are equally likely to prevail 
(this is the bistability condition) ; (ii) despite the interac- 
tion, individuals might decide to change their opinion in a 
random manner (stochastic terms in the dynamics); (hi) 
there is an external agent (advertising) that influences 
the opinion of the individuals (a forcing term). Given 
the aforementioned conditions, it should be of no sur- 
prise that a system size stochastic effect has been found 
in specific models. 

This result is quite general, and we now explain in de- 
tail a specific model which displays this effect. We have 
chosen a particularly simple opinion model developed by 
Kuperman and Zanette 66] . In this models each of the 
i = 1, . . . , N individuals of a society can have a binary 
opinion, say p,i — ±1, on a given topic. The individuals 
interact through a network of neighbors. We have con- 
sidered both the Watts- Strogatz small world networks 
as wel as the scale-free network, constructed with the 
Barabasi- Albert algorithm |5l| . 

The first ingredient of the dynamical evolution of the 
opinion is particularly simple and it is based upon a ma- 
jority rule: at a given time step t an individual i is se- 
lected, this individual then adopts the majority opinion 
held by his set of neighbors (there is no bias rule: in case 
of a tie, a random value for \ii is selected). This step, 
acting alone, would lead to a society holding a single 
opinion, either +1 or — 1, i.e, the dynamics is bistable. 
The selected value, +1 or —1, depends on the initial con- 
dition and the particular set of individuals chosen at each 
time step. 



The second ingredient reflects the existence of an exter- 
nal influence agent such as advertising: with probability 
ecos(wi)| the opinion of the randomly chosen individual 
i changes to sign(cos(wt)). Here, uj = 2n/T is the fre- 
quency with which the externally favored opinion changes 
with time. This periodic variation aims to reflect in some 
simplified manner the periodic variation of fashion. 

The third and final ingredient of the dynamics is ran- 
dom choice. Whatever the results of the two previous 
steps, individual i decides with a probability r\ to change 
his opinion to a randomly chosen value. 

After these three dynamical steps are taken consecu- 
tively (in the order specified above) by individual i, time 
increases by a time step t — ► t + 1/N. 

We look at the average opinion p = N^ 1 £\ /Ltj. In the 
absence of the periodic forcing, e = 0, it is found that p 
fluctuates stochastically between two values close to +1 
and — 1. In fact, the typical evolution, see the top panel 
of figure OH is very similar to the stochastic dynamics of 
a particle in a double well potential. A characteristic 
switching (Kramers) time can be defined as the aver- 
age time to jump from one favored opinion to the other. 
It is found that this time decreases with noise intensity 
r\ and increases exponentially with system size N. Al- 
though the dynamical rules do not allow us to identify a 
bistable potential function, it is clear that qualitatively 
the dynamics is such that a barrier between the two sta- 
ble opinions has to be overcome in order for the society 
to change opinion. 

This analysis confirms that we have at hand in this 
simple opinion model all the ingredients necessary for 
stochastic resonance: a bistable stochastic dynamics with 
a typical switching time Tk- As this time decreases with 
noise rate, it is possible to find a value of 77 such that half 
the external period coincides with the switching time, 
and hence the system responds optimally. This result 
was proven in reference 66] , although the first prediction 
for stochastic resonance in a model for opinion formation 
was obtained by Baninec in 68]. We show this result in 
figure ITU1 where we plot the spectral amplification factor, 
defined as R — 4e -2 | (e l 2lTt / T p(t)) | ,((...) denotes a time 
average). This is known to be a good measure of the 
response of the system to the forcing signal [6^. 

The fact that the switching time Tk depends on the 
system size N, makes it possible to achieve the resonance 
condition by varying N. This interesting result shows 
that the influence that the forcing signal has on the av- 
erage variable p depends on the size N. In other words, 
that the effect of advertising is maximum for a popu- 
lation of the right size, neither too large nor too small. 
This result is shown in figure lTTl where it is shown that the 
spectral amplification factor R has a maximum as a func- 
tion of the system size N. Although this result has been 
obtained within the framework of an specific model, it is 
important to realize that it wil also be present in other 
opinion models with the same basic ingredients. For in- 
stance, as shown in the figure lTTl the network structure is 
not a relevant variable here since the amplification factor 
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FIG. 9: We plot the time-evolution of the average opinion as a function of time (left column), and the power spectral density 
S (right column). The first row (a) shows the dynamics in absence of an external influence e = in the case N = 100. The 
last three rows, shows the results for different system sizes: N = 10, 100, 1000 (b, c, and d, respectively). The external signal 
(the sinusoidal thin line in the evolution plots) has a period T = 128. Note that this signal is better followed by a system 
with an intermediate size (panel c). This is also signalled by the largest peak of the spectral density at the driving frequency: 
fi = 27r/128. 



always displays a maximum at an optimal value of the 
system size. Of course, the exact value of the optimal 
size does depend on the network structure. 



IV. CONCLUSIONS 

In this paper, we have reviewed some situations in 
which system size can have a non-trivial role in some 
models of relevance to social systems. The first example 
was that of apparent phase transitions in the Galam's 
model of biased opinion. We showed in that example 
that there is a clear change of behavior in the sense that 
it is possible to identify a transition point separating the 
two possible final opinions. The transition point occurs 



at a pseudo-critical value p c for the initial fraction of sup- 
porters of the favored opinion. However, the value of p c 
depends on the system size in such a way that p c tends to 
zero when the system size N increases to infinity. There- 
fore, in a strict thermodynamic sense, the transition does 
no exists since it does not survive the thermodynamic 
limit. This is similar to what happens in one-dimensional 
models for the growth of rough surfaces, for which it can 
be proven that a change of behavior between a rough and 
a smooth phase occurs at some pseudo-critical value of 
the temperature. However, in accordance with some gen- 
eral theorems, there can not be a true phase transition 
in this one-dimensional system. 

Next we have unveiled a similar apparent phase tran- 
sition in Axelrod's model for culture formation. Here 
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FIG. 10: We plot the system response (measures through the 
spectral amplification lactor R) for different network topolo- 
gies in the Kuperman-Zanette model for opinion spreading. 
The different curves correspond to different network topolo- 
gies. The curves with open symbols correspond to a regular 
network (circles), a small world networks with rewiring proba- 
bility p — 0.20 (squares) and a fully random network, or p — 1 
(diamonds). The black triangles correspond to a Barabasi- 
Albert scale free network. The lines are only a guide for the 
eyes. The system parameters are N = 2 x 10 3 , T s = 128, 
e = 0.02. 
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FIG. 11: Spectral amplification factor as a function of N for 
different network topologies: gray symbols show the results 
for a small- world network with rewiring probability p — 0.2, 
black symbols are for a fully random network and white sym- 
bols correspond to the Barabasi- Albert network. In all the 
plots it is apparent an optimum response for an intermediate 
system size. The system parameters are as in the previous 
figures. 



the apparent phase transition corresponds to a change 
of behavior in the existence of global or multicultural 
states as a function of the noise intensity, as measured 
by the existence of cultural drift. Again, although a 
clear change of behavior can be observed, the transition 
does not exist in the true thermodynamic sense, since 
the pseudo-critical value of the noise intensity in regular 
lattices scales as r c oc l/Nh\N in two dimensions or as 
r c oc 1 /N 2 in one dimension. We have also considered the 
role of the structure of the social interactions by studying 
Axelrod's model in small world and scale free networks. 
It is found that in the Barabasi- Alert network, a pseudo- 
critical value q c for the heterogeneity parameter can be 
identified as separating the regimes of globalized and lo- 
calized states. In the thermodynamic limit, q c tends to 
zero and the transition disappears. 



Finally, we have presented a simple model for opin- 
ion formation based upon a majority rule in which it is 
shown that the effect that an external forcing (such as 
advertising) has on the society depends on the popula- 
tion size. A society of the proper size, neither too large 
nor too small, would react optimally to a weak signal. 



This work has stressed the non-trivial role that the 
system size has in the dynamics of social systems. There 
are changes of behavior which appear only in finite sys- 
tems. In a strict sense, these changes of behavior can not 
be called phase transitions since they disappear in the 
thermodynamic limit. However, one can not forget that 
social systems are never in the thermodynamic limit and 
that those effects can not be foreseen should one insist in 
taking this mathematical limit, as statistical physicists 
are routinely used to do. 
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